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Abstract 



An attempt to extract critical exponents 7, (3 and r from data on gold nuclei 
fragmentation due to interactions with nuclear emulsion at energies 4.0 A GeV and 
10.6 A GeV is presented. Based on analysis of Campi's 2nd charge moments, two 
subsets of data at each energy are selected from the inclusive data, corresponding to 



'liquid' and 'gas' phases. The extracted values of critical exponents for the selected 
data sets are in agreement with predictions of 'liquid-gas' model of phase transition. 
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1 Introduction 



Multifragmentation, a breakup of an excited nucleus into many intermediate 
mass fragments, has been discussed for almost twenty years in terms of sta- 
tistical mechanics; possible critical behavior was investigated. In its ground 
state nuclear matter behaves like a liquid. Mean field theory simulations [1] 
predict that nuclear equation of state resembles that of Van der Waals gas. 
Therefore existence of phase transition, spinodal instabilities and critical point 
are expected. However, nature of the possible phase transition is still under 
debate. Existing experimental data do not lead to a conclusive answer [2-5]. 

In many papers it was shown (e.g. [3]) that nuclear interactions undergo two 
stages. During the first stage, prompt nucleons are emitted from the colliding 
system and they carry out a large amount of available kinetic energy. They 
result from quasi-elastic and nonelastic collisions of projectile and target nu- 
cleons. Immediately after the collision, the remnant of the nucleus is in an 



excited state with temperature Tj. At the second stage, the excited remnant 
expands and cools evolving into neighborhood of the critical point on the 
temperature-density plane [6]. Then the hot and decompressed nucleon gas 
with temperature Tf condenses into many fragments. This last process is the 
multifragmentation. It was also shown that total charged fragment multiplic- 
ity is proportional to the temperature of the colliding system, both to 7* and 
Tf[6}. 

In the 90-ties there were many attempts to extract critical exponents in nuclear 
fragmentation from the experimental data, e.g. [7-10]. A method of charged 
moments invented by Campi [11] and supported by percolation theory [13] was 
commonly applied. These attempts did not take into account the fact that it is 
only the remnant of the nucleus that undergoes multifragmentation process, 
so all prompt particles (participating in the nucleus-nucleus collision) were 
included in the analysis. This approach raised a critique [14,15], that pointed 
out that: inclusion of prompt particles, assumption that fragment multiplicity 
is directly proportional to the temperature of the fragmenting system and 
using fragment charges Zi instead of fragment masses Ai are not justified and 
can significantly change the resulting values of critical exponents. 

Recently, the EOS collaboration [16] calculated critical exponents based on a 
high-statistics sample of fully reconstructed events from fragmentation of Au 
nuclei with energy 1.0 GeV/nucleon, taking into account the above criteria 
and using fragment masses, not charges. The obtained results are in agreement 
with the previously extracted values of the critical exponents [7]. Therefore 
the exclusion of prompt particles and performing analysis based on moments 
of mass distribution, instead of charge distribution, does not impact resulting 
values. This can be understood based on the fact that a large majority of 
prompt particles are fragments with charge Z — 1. Their impact on the second 
charge moments (that take square of the charges) is much smaller than that 
of heavier fragments. 

In this paper we present the analysis of data coming from interactions of pro- 
jectile gold nuclei of primary energies 4.0 and 10.6 GeV/nucleon with nuclear 
emulsion target. 



2 The experiment 

Stacks of BR-2 emulsion pellicles were irradiated with gold ion beam at the 
AGS accelerator at Brookhaven National Laboratory. The stacks were ori- 
ented so that the beam was parallel to the pellicles. Interactions were found 
during a microscope scanning along the primary tracks in order to obtain a 
sample with minimum detection bias. The two data sets consist of 448 events 
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at 4.0 GeV/nucleon and 884 events at 10.6 GeV/nucleon. In each event an- 
alyzed, multiplicities and emission angles of all produced particles and frag- 
ments of colliding nuclei were measured. In addition, charges of projectile nu- 
cleus fragments were determined. Singly-charged particles (released protons 
and produced pions) were distinguished unambiguously from heavier frag- 
ments. Charges of heavier fragments (Z > 2) of the projectile were measured 
using a photometric method with a CCD camera [17]. 

The experimental method of identification of prompt particles (i.e. nucleons 
directly participating in the collision) in emulsion experiments is not avail- 
able. We can however estimate the number of prompt particles as a difference 
between the total number of emitted protons N prot and the number of spec- 
tator protons N spec . N prot is determined using charge balance of the projectile 
fragments. We have estimated N spec as the number of singly charged particles 
that are emitted at the angle 6 < 28 , where 8q is the average proton emission 
angle given by relation 6 = 0.12/p, and p is the momentum of the nucleon 
prior to emission [18]. 

In this paper critical exponents were extracted using the Campi's charge mo- 
ments method. The charge moments in our analysis were calculated in two 
different ways. One of them is to find non-normalized charge moments taking 
into account all heavy [Z > 2) fragments, alpha particles and protons emitted 
from the vertex of interaction following Gilkes at al [7]. The other way is for 
normalized charge moments and excludes prompt protons following Elliott et 
al. [16]. The two ways were used to detect possible inconsistencies coming from 
the choice of the method. 



3 Charge moments 

Following Campi [11,12], we define the total charged fragment multiplicity m, 
as m = Nf + N a + N prot , where Nf denotes the number of fragments with 
charge Z > 3, N a is the number of emitted alpha particles and N prot is the 
number of emitted protons. 

The distance from the critical point for a given event can be properly measured 
by a difference between multiplicity m and the multiplicity at the critical point 
m c . The validity of this assumption rests on a linear dependence between 
temperature of the system T and the total multiplicity m, that was shown to 
be valid by Hauger at al. [3]. Therefore we introduce variable e: 

e = m — m c , (1) 
Fragments of a given charge, Zf, were counted on event-by-event basis to 
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determine fragment charge distributions Nz r The normalized charge distri- 
bution was defined as nz } = N Z} ./Za q , where Z^ is charge of the nucleus 
remnant undergoing fragmentation, given by a sum of charges of all fragments 
Zf > 2 plus total charge of spectator protons. 

In order to compare our results with theoretical calculations (the Fisher model), 
mass distributions should be used and should be normalized to the remnant 
mass A rather than to its charge Z^ . In emulsion experiments, fragment 
mass is not measurable; therefore it was assumed that fragment charge Zf is 
proportional to its mass Af. It was also assumed that on average, charge and 
mass distributions are equal N(Zf) = N(Af). In fact, the same assumption 
is made also by other authors, at least in specific ranges of charges. This is 
dictated by the difficulty in measuring masses of all heavier fragments emerg- 
ing from the interaction vertex. In case of electronic experiments masses of 
fragments only in some small mass range are measured. Therefore, we have to 
rely on the assumption of proportionality between mass and charge. 

Following Campi [11], we define the k-th moment of charge distribution as: 

M k (e) = J2n Zf (e)Z k f , (2) 

where the sum extends over all charged fragments except prompt protons in 
the 'gas phase', that is for e > 0. In the 'liquid phase',where e < 0, in addition 
the fragment with highest charge Z max is omitted from summation. The above 
procedure is motivated by the Fisher model. The bulk liquid of infinite volume 
is excluded from calculation on the 'liquid' phase. Similar procedure is carried 
out in percolation models, where percolating cluster of infinite size is excluded 
from calculation on the 'liquid' phase. 

With the above assumptions, and based on the Fisher model, in the thermo- 
dynamic limit the following relations are valid [16]: 

M 2 (e) ~ \e\~\ (3) 

n Zf (e) ~ ZJ T for e = (4) 

In addition, Bauer et al [19] have shown that: 

Z max (e) ~ t P for e < 0, (5) 

In these equations (3, 7 and r are the critical exponents. Z max (e) is the average 
charge of the largest fragment for a given value of the distance from critical 
point e. The above relations are valid in the neighborhood of the critical point. 
Far away from the critical point, the behavior of the system is dominated by 
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the mean field regime and these relations are not followed. Very close to the 
critical point, on the other hand, finite size effects come into play and M 2 
does not raise to infinity with e approaching (Equation 3), but achieves 
some maximum value. The critical exponents f3, 7 and r are not independent. 
A scaling relation between them exists: 



Practical calculations of the moments of charge distributions were based on 
averaging in the small bins of multiplicity m of event-by-event distributions 



where % is the index of an event, N denotes the total number of events in a 
given small range of e, M l k is a charge distribution moment for event i. 



4 Fluctuations close to the critical point 

One of the basic effects of the second order phase transition is the appear- 
ance of significant fluctuations in the neighborhood of the critical point, in a 
small range of temperatures T (or other parameter measuring distance from 
the critical point). Fluctuations grow as the critical point is approached and 
appear at increasingly large scales. In the case of a normal liquid, the effect 
is known as the critical opalescence: fluctuations of density of the liquid and 
sizes of gas bubbles result from vanishing latent heat of the phase transition as 
the critical point is approached. In the Fisher model of nuclear fragmentation, 
this is reflected by the divergence of the isothermal compressibility k t at the 
critical temperature T c : small variations of the pressure result in big density 
changes. In the neighborhood of the critical point, the volume and surface 
terms of the Gibbs free energy of fragment formation vanish and the fragment 
distribution is dominated by the power law of Equation 4 [20]. 

In case of multifragmentation, fluctuations can be analyzed using moments 
of charge distributions. Campi suggested a variable 72 linearly dependent on 
variance a 2 of the charge (Z) distribution: 



[16,7]: 
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Fig. 1. (72) as a function of total charged fragment multiplicity m at energies 4.0 
and 10.6 GeV/nucleon. The error bars show the dispersion of 72 in each bin. 

Figure 1 shows (72) as a function of multiplicity m for energies 4.0 and 10.6 
GeV/nucleon. As critical value of multiplicity m c is not identified at this stage, 
in all above calculations charge moments were calculated with Z max excluded. 
For multiplicities smaller than m ~ 30, (72) grows monotonically with mul- 
tiplicity, for multiplicities larger than m ~ 40, it falls down monotonically. 
Fluctuations of 72, reflected in large dispersion, are largest for multiplicities 
between 10 and 40. It is commonly assumed that strong maximum of (72) re- 
sults from undergoing a phase transition at critical point multiplicity m = m c . 
The region of multiplicities m < m c is called a 'liquid' phase, and events with 
m > m c are called to be in 'gas' phase. However, the presence of maximum 
of (72) is not a conclusive argument for appearance of phase transition. It 
was shown [16] that such maximum is observed also in systems that do not 
undergo a phase transition. Therefore Figure 1 will be treated as only a hint 
for possibility of phase transition in a specific range of multiplicities. It is also 
worth mentioning that 72 fluctuations at energies 4.0 and 10. 6A GeV agree 
well with each other. 



In order to closer analyze the fluctuations, the mean values of second moment 
of charge distribution (M 2 ) were plotted versus multiplicity in Figure 2. Similar 
characteristics to that of (72) is seen. Large dispersion of M 2 , especially in the 
'liquid' phase reflects strong fluctuations of experimental values. 

Figure 3 shows a scatter plot of second charge moments for individual events 
M\ versus multiplicity m at energy 10.6 GeV/nucleon. Strong fluctuations of 
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Fig. 2. Mean value of the second charge moment (M2) as a function of multiplicity 
m for energies 4.0 and 10.6 GeV/nucleon. The error bars show dispersion of M2 
values for each bin. 
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Fig. 3. Scatter plot of second charge moments for individual events M\ as function 
multiplicity m for energy 10.6 GeV/nucleon. 
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Fig. 4. Selection of 'liquid' and 'gas' groups of events in the M\ — m scatter plot. 



M\ for multiplicities 15 < m < 45 are clearly visible, which result in large 
dispersion of M 2 values (shown as error bars in Figure 2). For systems under- 
going a second order phase transition, it is expected that fluctuations rise in 
the neighborhood of the critical point. Therefore observation of large fluctu- 
ations quite far away from expected critical region should be attributed to a 
different physical phenomenon. For example, events with values of M\ > 4 and 
multiplicity m < 20 are fission-like (only 2 heavy fragments plus alphas and 
protons). Two distinct groups of experimental points are seen in Figure 3. One 
of these groups consists of events with multiplicity smaller than m = 35 and 
M2 < 0.9, the other — of events with multiplicities larger than m = 30 and 
largest M2 values for a given m. The former group, with small multiplicities 
and small M\ suggests that the fragment charge distribution must be domi- 
nated by one heavy fragment - i.e. the expected characteristics of the 'liquid' 
phase. On the other hand, the latter group of events, with large multiplic- 
ity and larger M2, suggests existence of a larger number of small fragments, 
thus it resembles the expected characteristics of the 'gas' phase. The events 
which do not belong to either of these groups may be thought of as nuclei in a 
'mixed' phase (see Section 7). The 'liquid' and 'gas' groups were selected from 
experimental data for further analysis. Figure 4 shows the selection criteria. 
The 'liquid' group of events consist of those with multiplicity smaller than 
m = 35 located below the solid line, while the 'gas' group consists of events 
with m > 30 and located above and to the right of the dashed line. Similar 
selection was performed for data with energy 4.0 GeV/nucleon - see Figure 5. 
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Fig. 5. Second charge moment M\ for individual events vs. multiplicity m for energy 
4.0 GeV/nucleon. Events selection criteria . 

Let us look at these two selected groups in order to check if they really can 
be treated as 'liquid' and 'gas' phases. For example, Figure 6 presents charge 
distributions as a function of heavy (Z > 2) fragment multiplicity, Nf, for the 
two selected groups of data. Fragment charge distributions are clearly differ- 
ent in both groups. The charge distribution for the group of low multiplicity 
events consists mainly of events with one heavy fragment emitted [Z > 45). 
If the second fragment is emitted, it has a small charge (Z < 3). Therefore, 
we will refer to this group of data as 'liquid' as for these events one large 
fragment remains after collision in analogy to the large drop of liquid. In most 
cases there are only alphas and protons accompanying the large heavy frag- 
ment. The second group of events has totally different charge distribution that 
consists of events with number of emitted heavy fragments ranging from 1 to 
9 and charges of fragments Z < 25. The majority of fragments in this group 
has charges Z < 15. We will call this group of events a 'gas' as it consists of 
events with many light fragments emitted. The analysis presented in subse- 
quent sections presents results on these selected sets of data unless explicitly 
stated otherwise. 



5 Critical exponent 7 



In order to extract the critical exponent 7 and the critical value of multiplicity 
m c the method known as 'gama matching' was used [19,16]. The outline of 
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Fig. 6. Charge distributions in events with different heavy fragment multiplicity Nf 
for the two selected groups of experimental data at energy 10.6 GeV/nucleon. 



this method is the following. A trial value of the critical multiplicity m c is 
chosen. For a given m c , a distribution of mean values of the moment (M 2 (e)) 
is determined as a function of distance from the critical point e = m — m c . 
Then the ranges in e are chosen to fit the power law (3) to the experimental 
data, separately for the 'gas' and 'liquid' phases. With fitting boundaries de- 
termined, the linear fit to the ln{M 2 (e)) versus ln\e\ is made to extract values 
of the slope 7 separately for 'gas' and 'liquid' phases. 

The fitting boundaries for the 'liquid' phase are easily determined. Figure 
7 presents a comparison of second moments of charge distribution for the 
selected 'liquid' phase and all available experimental data at m < m c . It is 
seen that the selection results in appearance of a very clear region of power 
law dependence. 

Similar comparison for the 'gas' phase is shown in Figure 8. The impact of 
selection of the 'gas' phase is not significant in this case. The region of the 
power law dependence is not clearly seen. This may be due to the inclusion of 
the charge of the largest fragment Z max in the 'gas' phase. The largest charge 
dominates the second moment, shifts it to the higher values and makes the 
power law less visible, probably due to the finite size effect. Thus, it may be 
interesting to check if the exclusion of Z max also in the 'gas' phase reveals a 
clear range of power law dependence of < M2 > versus e. The M2 moment 
for the 'gas' phase computed without Z max is shown in Figure 9. Open circles 
represent moments computed for 'gas' group of data, filled circles represent 
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Fig. 7. Mean value of second charge moment (M2) for a 'liquid' phase as a function 
of multiplicity m for energy 10.6 GeV/nucleon. Full squares show data for all events 
with m < m c , open squares represent data for selected 'liquid' events only. The line 
shows the fit used to determine critical exponent 7. 

those for all experimental data at m > m c , without selection. In both cases the 
region of power law behavior is now clearly visible. This region is much larger 
for selected 'gas' data compared to m > m c data, so it is now clear in which 
region of e the fit should be performed. In other words, Figure 9 was used 
as a guideline in determining fitting region in the 'gas' phase in the 'gamma 
matching procedure'. It is interesting that the value of 7 sas = 1.14 ± 0.05 
determined from moments of charge distribution with Z max excluded is very 
close to the value determined by the 'gamma matching procedure' (i.e. with 
Z max included). 

The above procedure was repeated for several trial values of the critical point 
m c . The value of the critical exponent 7 and critical point m c is found by 
demanding that \^g as — luquid\ takes a minimum value and ^ gas and ^u qu id agree 
with each other within statistical errors. Trial values of the critical point were 
selected in the range 25 < m c < 35. 

Results of the calculations for energy 10.6 GeV/nucleon are presented in Tables 
1 and 2. Table 1 shows results obtained based on normalized charge distribu- 
tion moments whereas Table 2 present results from calculation based on non 
normalized moments. It is easily seen that both procedures give exponents 
that agree with each other within statistical errors. Final value of exponent 
7 is calculated as a mean value of 7 9as and ^uquid for normalized moments. 
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Fig. 8. Mean value of second charge moment (M2) for a 'gas' phase as a function of 
multiplicity m for energy 10.6 GeV/nucleon. Full squares show data for all events 
with m > m c , open squares represent data for selected 'gas' events only. Plotted 
lines show fits to the 'gas' data (solid line) and m > m c data (dotted line). 
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Fig. 9. Mean value of the second charge moment (M2) for a 'gas' phase with ex- 
clusion of the largest charge Z max as a function of multiplicity m for energy 10.6 
GeV/nucleon. 
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Table 1 

Critical exponent 7 for 'gas' and 'liquid' phases for different choices of trial value of 
critical point m c for energy 10.6 GeV/nucleon used in gamma matching procedure. 
Normalized moments 



trial m c 


^/liquid 


'Jgas 


Illiquid 1gas\ 


26 


0.80 ± 0.00 


1 


45 ± 0.06 


0.63± 0.10 


28 


0.96 ±0.11 


1 


32 ± 0.05 


0.35 ±0.12 


29 


1.04 ±0.11 


1 


26 ± 0.05 


0.22 ±0.12 


30 


1.12 ±0.12 


1 


20 ± 0.05 


0.08 ±0.13 


31 


1.20 ±0.13 


1 


14 ±0.05 


0.05 ±0.14 


32 


1.27 ±0.14 


1 


08 ± 0.05 


0.19 ±0.15 


34 


1.42 ±0.15 





96 ± 0.04 


0.46 ±0.16 



Table 2 

Critical exponent 7 for 'gas' and 'liquid' phases for different choices of trial value of 
critical point m c for energy 10.6 GeV/nucleon used in gamma matching procedure. 
Non-normalized moments. 



trial m c 


I liquid 


1 gas 


lOgcts ^liquid] 


26 


0.78 ± 0.07 


1.30 ±0 


11 


0.48± 0.13 


28 


0.93 ± 0.09 


1.22 ±0 


10 


0.20 ±0.14 


29 


1.02 ±0.10 


1.17±0 


10 


0.15 ±0.14 


30 


1.09 ±0.10 


1.13 ±0 


10 


0.03 ±0.14 


31 


1.17 ±0.10 


1.08 ±0 


10 


0.09 ±0.14 


32 


1.24 ±0.12 


1.03 ±0 


09 


0.21 ±0.14 


34 


1.39 ±0.13 


0.94 ± 


08 


0.45 ±0.15 



The critical point is determined at m c = 31 ± 2 and the critical exponent 
7 = 1.17 ±0.09. 

The analogous results of the 'gama matching' procedure for energy 4.0 GeV/nucleon 
are given in Tables 3 and 4. The determined critical point is m c = 30 ± 2 and 
the critical exponent 7 = 1.15 ± 0.09. Figure 10 shows values of \^u qu id — l gas \ 
plotted as a function of trial values of critical multiplicity at energy 10.6 
Gev/nucleon. Figure 11 presents results of fitting procedure for / yii qu id and 
7 9as for critical multiplicity value m c = 30 for the same energy. 
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Fig. 10. Difference 
10.6 GeV/nucleon. 



Illiquid ~~ lgas\ as a function of trial multiplicity m c for energy 
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Fig. 11. Mean value of second charge moment (M2) for a 'liquid' and 'gas' as a 
function of multiplicity m for energy 10.6 GeV/nucleon. The plotted lines show fits 
to experimental data used to determine critical exponent 7. 
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Table 3 

Critical exponent 7 for 'gas' and 'liquid' phases for different choices of trial value of 
critical point m c for energy 4.0 GeV/nucleon used in gamma matching procedure. 
Normalized moments. 



trial m c 


^/liquid 


'Jgas 


I "Igas ^liquid \ 


26 


0.83 ± 0.09 


1.32 ±0.09 


0.49 ±0.13 


28 


1.01 ±0.11 


1.22 ±0.09 


0.21 ±0.14 


29 


1.10 ±0.12 


1.17 ±0.08 


0.08 ±0.14 


30 


1.18 ±0.13 


1.12 ±0.08 


0.06 ±0.14 


31 


1.27 ±0.14 


1.07 ±0.08 


0.20 ±0.15 


32 


1.35 ±0.15 


1.01 ±0.07 


0.33 ±0.16 


34 


1.52 ±0.16 


0.91 ±0.07 


0.61 ±0.17 



Table 4 

Critical exponent 7 for 'gas' and 'liquid' phases for different choices of trial value of 
critical point m c for energy 4.0 GeV/nucleon used in gamma matching procedure. 
Non-normalized moments. 



trial m c 


^liquid 


"fgas 


I "Igas ^liquid I 


26 


0.77 ±0.07 


1.53 ±0.08 


0.75 ±0.11 


28 


0.96 ± 0.08 


1.42 ±0.07 


0.46 ±0.11 


29 


1.05 ±0.09 


1.36 ±0.07 


0.31 ±0.11 


30 


1.13 ±0.09 


1.30 ±0.07 


0.17 ±0.11 


31 


1.22 ±0.09 


1.24 ±0.07 


0.02 ±0.11 


32 


1.30 ±0.10 


1.18 ±0.06 


0.13 ±0.12 


34 


1.47 ±0.11 


1.06 ±0.06 


0.41 ±0.14 



6 Critical exponents (3 and r 



After determination of the critical multiplicity m c and ranges of e in which 
to fit the critical exponents, determination of the exponents (3 i r is straight- 
forward. Based on Equation 5, the mean value of \nZ max was plotted as a 
function of In \m — m c \ (Figure 12). The slope of the linear fit to the plotted 
relation gives the critical exponent f3. A linear fit was made within the fitting 
boundaries determined during the 'gamma matching procedure'. For energy 
10.6 GeV/nucleon the determined value of (3 = 0.33±0.01 with x 2 / n df = 1.66, 
and for 4.0 GeV/nucleon f3 = 0.34 ± 0.01 with X 2 /ndf = 0.99. 
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Fig. 12. Mean Value of ln(Z max ) as a function of distance from critical point e for 
'liquid' phase 10.6 GeV/nucleon. 



The critical exponent r was determined from equation [7,16]: 



Aln(M 3 ) _ r-4 
Aln(M 2 ) ~ r-3' 



(9) 



Figure 13 presents ln(M 3 ) versus ln(M 2 ) for the 'gas' phase at energy 10.6 
GeV/nucleon together with the linear fit. The slope of this linear fit is used 
to determine the r exponent, which is r = 2.11 ± 0.05 for 10.6 GeV/nucleon 
and r = 2.12 ±0.04 for energy 4.0 GeV/nucleon. To verify the consistency, the 
exponent r also was determined from Equation 4. This Equation is supposed 
to be valid at the critical point, but in order to have sufficient statistics, data 
with 25 < m < 35 were used. Normalized charge multiplicity distribution was 
calculated and plotted in Figure 14. The exponent r is given by the slope of 
the linear fit to the data points for charges from Z = 6 to Z = 16. For charges 
smaller than 6, the assumptions of the Fisher model are not valid [22] and for 
charges larger than 16 experimental data statistics is too small. The resulting 
value r = 2.19 ± 0.33, with reduced x 2 — 1-2 is in agreement with previously 
determined value. 
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Fig. 13. Relation between charge moments ln(M3) and ln(Ms) for gas phase 10.6 
GeV/nucleon. 
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Fig. 14. Normalized charge distribution for event with multiplicity m within a range 
25 < m < 35. 
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Table 5 



Critical exponents 7, f3, r for energy 4.0 GeV/nucleon and 10.6GeV/nucleon 



Data 


7 


(3 


r 


Au-Em 4. OA GeV norm 


1.15 ±0.09 


0.34 ±0.01 


2.12 ±0.04 


Au-Em 10.6A GeV norm 


1.17 ±0.09 


0.33 ±0.01 


2.11 ±0.05 


Au-Em 4.0A GeV 


1.23 ±0.09 


0.36 ± 0.02 


2.15 ±0.04 


Au-Em 10.6A GeV 


1.11 ±0.09 


0.33 ± 0.02 


2.16 ±0.05 


Au-C 1.0A GeV (EOS [16]) 


1.40 ±0.10 


0.29 ± 0.02 


2.14 ±0.06 


Au-Em 10.6A GeV (EMU01 [10]) 


0.86 ± 0.05 


0.25 ±0.02 


2.23 ±0.05 


Au-Em 10.6A GeV (KLMM [8]) 




0.19 ±0.02 


1.88 ± 0.06 


Liquid-gas model 


1.23 


0.33 


2.21 


Percolation model 


1.80 


0.41 


2.18 



7 Discussion 

Table 5 summarizes the determined values of critical exponents 7, (3 and r at 
energies 4.0 and 10.6 GeV/nucleon both for normalized and non-normalized 
charge moments. For comparison, exponents obtained by the EOS, EMU01 
and KLMM experiments are also shown, as well as the values predicted by 
the liquid-gas and percolation models. The data on the (5 and 7 exponents from 
Table 5 are plotted in Figure 15. The values of (5 and 7 exponents determined in 
this work are very close to those expected for the 'liquid-gas' phase transition, 
both for 4.0 and 10.6 GeV/nucleon data. We note that the critical exponents 
discussed here were determined from the selected 'liquid' and 'gas' groups of 
events. Earlier analyses, without such a selection, gave values of the exponents 
which do not agree with any of the models. 

Values of the critical exponent r predicted by the percolation and liquid-gas 
models are not too different. As experimental errors are large compared to the 
difference between the predicted values of r, this exponent cannot be used to 
discriminate models of multifragmentation. 

Presented in this work a suggestion that the Fisher's liquid drop model prop- 
erly describes multifragmentation process, is in agreement with recently pub- 
lished results of the ISiS collaboration [23] , confirming that the Fisher's scaling 
law is followed by experimental data. In addition, Mader et al. [24] have shown 
a similarity between predictions of the Fisher model and clusterisation process 
in the three-dimensional Ising model. Therefore, multifragmentation can be 
interpreted in a similar way as condensation of liquid drops in equilibrium 
with the bulk liquid. 
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Fig. 15. Comparison of determined and predicted values of critical exponent 7 and 

(3. 



One can try to estimate a critical temperature of the system by using a relation 
between multiplicity and temperature of the system given in [3,6]. Estimation 
of the initial temperature of the system after collision, Tj, was based on Fermi 
gas model and does not take into account the expansion of the system. The 
final temperature of the system at breakup, Tf, was estimated based on a re- 
lation between multiplicity m and temperature T for isotopic-yield-ratio ther- 
mometer given in [3]. Our values of critical multiplicity (m c = 30 and m c = 31 
at 4.0 and 10.6 GeV/nucleon, respectively) correspond to T = 9.3 ± 0.7 MeV 
and Tf = 5.4 ± 0.2 MeV. The error of T reflects uncertainty of the input pa- 
rameters in the Fermi gas model. These T and Tf values bracket the critical 
temperature of the system. This estimated range for the critical temperature 
agrees with predictions of the Fisher model if we take into account the scaling 
of critical temperature with the finite system size [16]. The estimation of the 
critical temperature relies on the assumption that relation between multiplic- 
ity m and temperature of the system established for energy 1.0 A GeV applies 
also for higher energies 4.0 and 10.6 GeV/nucleon. 

The above results, together with results of ISiS collaboration [23] suggest that 
not all multifragmentation events undergo a second order phase transition. 
The data excluded from our analysis presented above may be interpreted as 
events in which multifragmentation occurs far away from the critical point. 
These events could undergo a first order phase transition so that 'gas' and 
'liquid' coexist inside the nucleus as suggested by experimental analysis of 
GSI data [2]. 
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As noted earlier, the excited nucleus after a collision evolves into the neigh- 
borhood of the critical point on the density-temperature plane. If multifrag- 
mentation occurs close to the critical point, critical exponents characterizing 
second order phase transition occurring in the neighborhood of the critical 
point should be observed. If for some events multifragmenation occurs far 
away from the critical point and they are included in the analysis, the ob- 
served values of critical exponents could be altered so that the whole picture 
of the physical process is obscured. In such situations proper selection of events 
should help obtain true values of critical exponents. 

Events excluded from our analysis show very strong fluctuations of the sec- 
ond moment of the charge distribution significantly far away from the critical 
point, especially in the liquid phase. Large fluctuations could result from coex- 
istence of two different phases with different properties inside a single nucleus. 
Below critical multiplicity m c (proportional to temperature of the system) two 
separate phases with significantly different properties could coexist in nucleus 
as is suggested by mean field [25,26,1] or canonical model calculations [27]. 
The coexistence of the two phases may take place in wide ranges of temper- 
ature and pressure. This could explain large fluctuations of second moment 
of charge distribution far away from the critical point. It is also important 
to stress that moment M 2 is proportional to the isothermal compressibility of 
the system. Therefore, large fluctuations of M 2 correspond to fluctuations of 
the compressibility kt and to fluctuation of the density of the system. Large 
fluctuations of M 2 observed in 'liquid phase' (Figure 3) may be interpreted as 
resulting from large density difference between the two coexisting phases. 

In summary, an attempt to extract critical exponents 7, (3 and r was per- 
formed, using data coming from interactions of gold nuclei with nuclear emul- 
sion at energies 4.0 A GeV and 10.6 A GeV. To extract the exponents, two 
subsets of data with characteristics similar to that of 'gas' and 'liquid' phases 
were selected, based on analysis of Campi's 2nd charge moments. The ex- 
tracted values of the critical exponents for the selected data sets are in agree- 
ment with predictions of liquid-gas model of phase transition. The same analy- 
sis performed without the selection of 'gas' and 'liquid' samples favors neither 
percolation nor liquid-gas model of phase transition. A suggestion is made 
that data excluded from the above mentioned samples represents events where 
phase transition, if any, occurs far from the critical point. 
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